Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 20 January 2013 



(MN style file v2.2) 



Dynamics of massive stellar black holes in young star 
clusters and the displacement of ultra-luminous X-ray 



sources 



M. Mapelli\ E. Ripamonti\ L. Zampieri^, M. Colpi-*^ 

^ Universitd Milano Bicocca, Dipartimento di Fisica G.Occhialini, Piazza della Scienza 3, 1-20126, Milano, Italy: inapelliOmib.infii.it 
^ INAF-Osservatorio astronomico di Padova, Vicolo deWOsservatorio 5, 1-35122, Padova, Italy 



20 January 2013 



o ; 



O: 
6 



00 

p 
o 



> ■ 1 INTRODUCTION 



X 



ABSTRACT 

In low-metallicity environments, massive stars might avoid supernova explosion and 
directly collapse, forming massive 25 — 80 Mq) stellar black holes (MSBHs), at the 
end of their life. MSBHs, when hosted in young massive clusters, are expected to form 
binaries and to strongly interact with stars, mainly via three-body encounters. We 
simulate various realizations of young star clusters hosting MSBHs in hard binaries 
with massive stars. We show that a large fraction 44 per cent) of MSBH binaries 
are ejected on a short timescale (< 10 Myr). The offset of the ejected MSBHs with 
respect to the parent cluster is consistent with observations of X-ray binaries and ultra- 
luminous X-ray sources. Furthermore, three-body encounters change the properties of 
MSBH binaries: the semi- major axis changes by < 50 per cent and the eccentricity of 
the system generally increases. We shortly discuss the implications of our simulations 
on the formation of high-mass X-ray binaries hosting MSBHs. 

Key words: black hole physics - stars: binaries: general - galaxies: star clusters: 
general ~ X-rays: binaries 



According to models of stellar evolution (Fryer 1999; Heger 
et al. 2002, 2003), a star with a final mass^ rrifln higher than 
mth ~ 40 M0 can avoid supernova explosion and directly 
collapse into a black hole (BH). In this case, it is reasonable 
to expect that the mass of the remnant is comparable to 
the final mass of the progenitor star (or at least more than 
half of it). Therefore, BHs may form, via this channel, with 
mass higher than in the standard scenario of supernova ex- 
plosion and fallback. According to Belczynski et al. (2010), 
the masses of BHs formed via direct collapse may be as high 
as ~ 80 Mq. In the following, we will refer to massive stel- 
lar BHs (MSBHs) to indicate BHs with mass 25 - 80 M© 
formed via direct collapse. 

The final mass of a star likely depends on metallicity, as 
mass losses by stellar winds correlate with metallicity (see 
Bresolin & Kudritzki 2004 for a review). Massive stars with 
(about) solar metallicity cannot have a final mass as high as 
~ 40 Mq , because they lose a large fraction of their mass by 
stellar winds. Thus, solar- metallicity stars never form BHs 



^ We name 'final mass', ma^, of a star the mass bound to the 
star immediately before the collapse. 



with mass higher than ~ 10 — 20 M0. Instead, stars with suf- 
ficiently low metallicity might have final masses ^ 40 M0 
and directly collapse into MSBHs. The metallicity thresh- 
old for the formation of MSBHs depends on the adopted 
model of stellar wind and stellar evolution. For example, as- 
suming the model by Belczynski et al. (2010), based on the 
stellar- wind model by Vink, de Koter & Lamers (2001), the 
threshold metallicity is Zth ~ 0.4 Zq. We stress that the de- 
pendence of stellar winds on metallicity is quite uncertain, 
especially for massive stars, as it relies upon observations 
of a few massive stars, and it is drawn from models gener- 
ally neglecting binary evolution and rotation (e.g., Bresolin 
& Kudritzki 2004, and references therein). For this reason, 
there is an uncertainty of ^ 0.2 dex on Zth- 

Recent papers (Mapelli, Colpi & Zampieri 2009; 
Zampieri & Roberts 2009; Mapelh et al. 2010a, 2010b) 
indicate that MSBHs born from the direct collapse of 
massive low-metallicity stars might account for a fraction 
of ultra-luminous X-ray sources (ULXs). ULXs are non- 
nuclear point-like sources with (isotropic) X-ray luminos- 
ity Lx ^ 10^® erg s"^ (see Mushotzky 2004 for a review, 
and references therein), corresponding to the Eddington lu- 
minosity of a ^7 Mq BH. The mechanism that powers 
the ULXs is still unknown, although various scenarios have 
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been proposed. Apart from the aforementioned interpre- 
tation (i.e., that a fraction of ULXs may be powered by 
MSBHs formed from the evolution of massive metal-poor 
stars), ULXs could be associated with high-mass X-ray bi- 
naries (HMXBs) powered by stellar-mass black holes (BHs) 
with anisotropic X-ray emission (e.g. King et al. 2001) or 
with super-Eddington accretion rate/luminosity (e.g. Begel- 
man 2002; King & Pounds 2003; Socrates & Davis 2006; 
Poutanen et al. 2007), or emitting via a combination of 
the two mechanisms (e.g. King 2008). Recently, Linden 
et al. (2010, hereafter LIO) suggested that a fraction of 
ULXs may be explained with Roche-lobe overflow (RLO) 
HMXBs in moderately low-mctallicity environments, under- 
going mildly super-Eddington accretion. In their scenario, 
the bulk of HMXBs hosts 10-15 M© BHs, but up to ~ 20 
per cent of HMXBs hosts BHs with mass > 20 Mq . In ad- 
dition, ULXs could also be associated with HMXBs pow- 
ered by intermediate-mass BHs (IMBHs), i.e. BHs with mass 
100 Mq < ruBH < 10^ Mq (see van der Marcl 2004 for a re- 
view). Among the aforementioned hypotheses, beamed emis- 
sion is problematic for those ULXs surrounded by isotropi- 
cally ionized nebulae (PakuU & Mirioni 2002; Kaaret, Ward 
& Zezas 2004a) and for the brightest ULXs. Similarly, super- 
Eddington emission requires uncommonly high violations of 
the Eddington limit for the brightest ULXs. On the other 
hand, the very existence of IMBHs lacks observational ev- 
idences and IMBHs are not needed to explain the obser- 
vational properties of most of the ULXs (e.g. Gongalves & 
Soria 2006; Stobbart, Roberts & Wilms 2006; Copperwheat 
et al. 2007; Roberts 2007; Zampieri & Roberts 2009). Fur- 
thermore, ULXs might actually be a mixed bag of different 
objects, including sources of different nature. 

In the following, we will focus on the scenario proposed 
by MapelU et al. (2009) and by Zampieri & Roberts (2009), 
according to which a fraction of ULXs might be powered 
by MSBHs, accreting in HMXBs. This idea is supported 
by the fact that most ULXs arc located in galaxies with 
a high star formation rate (SFR, e.g. Irwin, Bregman & 
Athey 2004) and are often associated with low-metallicity 
environments (Pakull & Mirioni 2002; Zampieri et al. 2004; 
Soria et al. 2005; Swartz, Soria & Tennant 2008; but sec 
Winter, Mushotzky & Reynolds 2007 for a different result'^). 
Considering a sample of 66 late-type galaxies, Mapelli et al. 
(2010a,b) showed that the number of observed ULXs per 
galaxy (iVuLx) scales almost linearly with the SFR and that 
-^ULx/SFR (i.e. the number of observed ULXs per galaxy, 
normalized to the SFR) anti-correlates with the metallicity. 
Finally, the expected number of MSBHs per galaxy (A^bh), 
derived in Mapelli et al. (2010a,b), correlates with A'ulx- 

Thus, the idea that MSBHs are the engine of most 
ULXs agrees, under many respects, with the data. However, 
many issues need to be addressed, in order to check this sce- 

2 Winter et al. (2007) derive the metallicity of 16 ULXs from 
XMM Newton spectra and find approximately solar abundances 
for all of them. However, X-ray photo-ionization edges are known 
to be highly uncertain metallicity indicators. The abundances de- 
rived by Winter et al. (2007) are much higher (0.2 — 1 dex) than 
those derived, for the same galaxies, from spectra of HII regions 
using the P-calibration method (see, e.g., Pilyugin, Vilchez & 
Contini 2004). For a more detailed discussion about metallicity 
estimators, see, e.g., Kennicutt, Bresolin & Garnett (2003). 



nario. First, observations indicate that ULXs are generally 
found close to star forming regions and to young star clus- 
ters (YSCs), but often displaced from them, by a distance 
« 0.1 — 1 kpc (e.g. Zezas et al. 2002; Swartz, Tennant & 
Soria 2009; Berghca 2009; Swartz 2010). The same seems to 
occur for bright (> 10^^ erg s~^) X-ray sources in starburst 
galaxies (e.g. Kaaret et al. 2004b). This fact has been gener- 
ally interpreted as the indication that bright X-ray sources 
are powered by runaway binaries, i.e. by binaries which were 
ejected from the parent cluster because of a natal kick (e.g. 
Sepinsky, Kalogera & Belczynski 2005; Zuo & Li 2010) or a 
close encounter (e.g. Kaaret et al. 2004b; Berghea 2009). The 
amount and the distribution of natal kicks arc very uncertain 
for BHs (e.g., Brandt & Podsiadlowski 1995; Brandt, Podsi- 
adlowski & Sigurdsson 1995; Gualandris et al. 2005; Fragos 
et al. 2009, and references therein). It is reasonable to as- 
sume that natal kicks in the absence of supernova are smaller 
than in the supernova scenario (e.g. Fryer & Kalogera 2001; 
LIO). If the natal kick is lower than the escape velocity from 
the parent cluster, it may not be sufficient to eject the MSBH 
from the parent cluster (e.g., LIO). LIO even assume that no 
natal kick is present in the direct collapse scenario, because 
of the absence of supernova explosion. While the role of na- 
tal kicks needs to be investigated in more detail, we focus 
here on the effects of recoil velocities induced by three-body 
encounters: are they able to eject MSBHs? 

A second issue is related to the evolutionary path for the 
formation of HMXBs with MSBHs. LIO suggest that a kick 
is almost necessary for a binary to enter the RLO-HMXB 
regime. According to their model, MSBHs can hardly enter 
the RLO-HMXB regime, as they have no natal kicks. How- 
ever, the model by LIO does not include three-body encoun- 
ters. In this paper, we want to study the effects of three-body 
encounters on the orbital parameters of a MSBH binary and 
on the possibility that such binary enters the RLO-HMXB 
regime. 

Therefore, the aim of this paper is to check whether 
there may be an alternative source of kick, different from 
the natal one, for MSBH binaries. In particular, we check, 
via direct N-Body simulations, the occurrence and the 
consequences of kicks by three-body interactions involving 
MSBHs in YSCs. 

We do not consider only MSBHs, but, for comparison, 
also intermediate-mass BHs (IMBHs). The plan of the paper 
is the following. In Section 2, we describe the simulations. 
In Section 3, we discuss the evolution of the radial distribu- 
tion of the simulated MSBHs. In Section 4, we mention the 
effects of the kick on the orbital parameters of the involved 
binaries. In Section 5, we present simulations of IMBHs em- 
bedded in YSCs, for comparison with MSBHs. In Section 6, 
we compare our simulations with the available data in the 
literature. In Section 7, we report our conclusions. 



2 SIMULATIONS 

The simulations were done using the Starlab^ software en- 
vironment (see Portegies Zwart et al. 2001), which allows to 



^ http://www.sns.ias.edu/~ starlab/ 
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integrate the dynamical evolution of a star cluster, resolving 
binaries and three-body encounters. 

We make 25 realizations of isolated star clusters having 
a King profile with adimensional central potential Wo =5 
(King 1966). Each cluster hosts Nstar ~ 5000 stars, corre- 
sponding to a total mass « 3 — 4 x 10'^ Mq. We set an 
initial binary fraction /bin ~ 0.1. Stars are generated ac- 
cording to the Salpeter (1955) initial mass function (IMF), 
in a mass range from 0.08 Mq up to 120 Mq. The virial 
radius of each cluster is ~ 1 pc. Thus, the clusters have 
a half-mass relaxation time trix ~ 10 Myr and a core col- 
lapse time scale tec — 0.2 t^ix ~ 1 Myr (see Portegies Zwart 
& McMillan 2002). We choose to simulate clusters with the 
same global properties to reduce the considered parameter 
space. No external tidal field is added. Stellar and binary 
evolution were switched off during these runs, to avoid in- 
troducing too many parameters in this first study. The ef- 
fects of binary evolution will be taken into account in paper 
II. Fig. 1 shows the initial one-dimensional profile of one of 
the realizations, compared with the profile at f = 10 Myr. 
According to Mapelli et al. (2010a), clusters with these prop- 
erties are expected to host f» 1 MSBH, if their metallicity is 
~ 0.1 Zq. Therefore, we generate one MSBH in each simu- 
lated cluster. The mass of the MSBH is fixed to be 50 Mq, 
as fiducial value. The MSBH is always generated in a binary 
system, with a massive companion. The assumption that 
the MSBH is always in a binary is quite realistic, as the pro- 
genitor of the MSBH and the MSBH itself are among the 
most massive objects in the parent cluster (see, e.g., Kulka- 
rni. Hut & McMillan 1993). The mass of the companion 
is uniformly distributed (following the standard settings of 
Starlab) between a minimum mass miow and a maximum 
mass mhigh. We require that miow — 10 M© and mhigh = 50 
Mq. In fact, theoretical studies (e.g., Portegies Zwart, Dewi 
& Maccarone 2004; Patruno et al. 2005; Madhusudhan et al. 
2006) indicate that, for masses of the donor m2 ^ 10 Mq, 
luminosities in the ULX range cannot be attained during 
main sequence (MS), except for very short outbursts a 
few days every few months of quiescence). We do not con- 
sider ULXs with lower-mass evolved donors, because we are 
interested only in young systems {t ^ 10 Myr), where stars 
less massive than ~ WMq did not yet evolve off the MS 
and transfer mass steadily on the MSBH. Furthermore, we 
note that if the mass of the donor is larger than the mass 
of the BH, the system will have an unstable mass transfer 
(Rappaport, Podsiadlowski & Pfahl 2005). For this reason, 
we avoid companion masses larger than the MSBH mass. 
The initial position and velocity of the centre of mass of the 
MSBH binary are drawn from the same distributions as the 
other stars in the cluster (see Portegies Zwart et al. 2001). 

The initial semi-major axis and the initial eccentricity 
of the MSBH binary follow the same distributions adopted 
for 'normal' binaries (i.e. a thermal distribution for eccen- 
tricities and a logarithmic distribution for semi-major axes, 
Portegies Zwart et al. 2001). We then exclude all MSBH 
systems with semi-major axis a > 10 A.U. and a < 0.1 
A.U. and all systems that were already ionized or that ex- 
changed companion after one timescale (corresponding to 
~ 0.25 Myr). We consider only systems with a > 0.1 A.U., 
to exclude binaries whose semi-major axis is smaller than 
the tidal radius rtid of the companion (for the binaries in 
our sample, rtid = 0.03 — 0.1 A.U.). Furthermore, we select 




Figure 1. Surface density of the cluster (S), normalized to 
the central surface density value So, as a function of the two- 
dimensional projected distance r. rs^ale is the scale radius and is 
~ 1 pc for all simulated clusters. Open circles (blue on the web): 
initial conditions; filled squares (red on the web): evolution after 
10 Myr. 



only binaries with a < 10 A.U., because we are interested 
only in systems that are sufficiently hard to survive for most 
of the cluster's life''. Therefore, our sample is not represen- 
tative of all possible binaries hosting MSBHs, but only of 
hard systems. We will consider softer binaries, whose dy- 
namics is completely different, in a forthcoming paper. We 
also stress that our simulations start when the MSBH is al- 
ready formed and do not make any assumptions about the 
dynamical and stellar evolution before the formation of the 
MSBH. Therefore, the simulated hard MSBH binaries are 
not necessarily primordial, but can have a dynamical ori- 
gin (via exchanges) . Furthermore, our selection of the initial 
semi-major axes provides initial periods that are consistent 
with the expected distribution of periods of hard BH binaries 
in dense star clusters (e.g., figure 2 of Belczynski et al. 2004). 
Each cluster realization has been integrated for 10 Myr, i.e. 
approximately the lifetime of a 15 Mq star. This guarantees 
that the companion of the MSBH is still unevolved in most 
cases. 

For comparison with the runs including MSBHs, we also 
simulated a set of 20 clusters with the same properties as 
before but hosting an IMBH (with mass 300 Mq) instead 
of a MSBH. The IMBH is always generated in a binary, 
with similar constraints to those adopted for MSBHs (i.e. 
companion mass between 10 and 120 Mq and semi-major 
axis such that the binary is hard). 



* A binary is hard when a < G mi m2 (mi -I- m2 + 
rti'i)/ [(mi -I- m2) ms cr^j , where G is the gravitational constant, 
mi and m2 are the mass of the primary and of the secondary, 
respectively, ma is the average mass of a cluster star, a is the 
velocity dispersion of the cluster. For our systems, the MSBH 
binary is hard if a ^ 10 A.U. 
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Table 1. End states of the simulated MSBHs that are ejected 
before t = 10 Myr. 



State 


Percentage^ 


Ionization 


0.0 


Binary 


81.8 


Triple 


18.2 


Exchange 


27.3 



^The percentage refers to the total of ejected MSBHs before 
J = 10 Myr (i.e. 11 cases). 
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3 CLOSE ENCOUNTERS AND EJECTION 

Fig. 2 shows the radial projected distribution of MSBHs at 
t = (empty black histogram) and at t = 10 Myr (filled 
histogram, red on the web). The position of each MSBH in 
Fig. 2 comes from a different run. We assume that a MSBH 
is definitely ejected from its parent cluster, when its pro- 
jected distance from the centre of the cluster is larger than 
4rscaio (i.e. ~ 4 times the virial radius, which corresponds 
approximately to the tidal radius of the clusters) and when 
its three-dimensional velocity is larger than 5 km s~^ (i.e., 
the escape velocity from the parent cluster). According to 
this criterion, 11 MSBHs are ejected at t < 10 Myr, i.e. 44 
per cent of the simulated MSBHs. Before the ejection, the 
MSBHs had a strong interaction with a single star, whose 
mass ranges from ~ 24 M© to ~ 61 M©. The ejected MSBHs 
are still in binaries (82 per cent, i.e. 9 runs) or even in triple 
systems (18 per cent, i.e. two runs), as it is reported in Ta- 
ble 1. This is a consequence of the fact that we consider 
only hard binaries. Furthermore, a non-negligible fraction 
of ejected MSBHs (27 per cent) exchanged their companion 
with a different star before the ejection. Table 2 shows the 
analogous statistics for MSBHs that, at t = 10 Myr, are still 
bound to the cluster. The percentage of ionizations among 
MSBHs that remain in the cluster is still very low, but the 
percentage of (dynamically unstable) triple systems is much 
larger than that of binaries. This means that most of MSBHs 
that remain inside the cluster are still interacting. Interest- 
ingly, no exchanges occur among the MSBHs that remain 
inside the cluster. 

Let us now focus on the properties of ejected MSBH sys- 
tems. Fig. 3 shows the distribution of recoil velocities (uroc) 
of the ejected MSBHs. v^cc spans from 5 km s"'^ up to 120 km 
s~^: it is, on average, smaller than the natal kick velocities 
assumed for stellar-mass BHs in LIO (see e.g. their fig. 6), 
but it is a viable substitute to natal kick for BHs born from 
direct collapse. Furthermore, velocities larger than ~ 100 
km s~^ are likely to disrupt binaries (e.g., Berghea 2009). 

Fig. 4 shows the distribution of times of MSBH ejection 
from the parent cluster. We calculate it with an accuracy 
of ~ 0.25 Myr. Most of ejections occur between 2.5 and 6.5 
Myr, when the massive companion of the MSBH is expected 
to be on the MS. 



Figure 2. Empty black histogram: two-dimensional radial pro- 
jected distribution of MSBHs at time t = 0. Filled histogram 
(red on the web): two-dimensional radial projected distribution 
of MSBHs at time t = 10 Myr. The empty black histogram and 
the filled histogram refer to left-hand j/— axis (i.e. the number 
of simulated MSBHs, A'^bh) and to the bottom x'— axis (i.e. the 
two-dimensional distance r of the MSBH from the centre of the 
parent cluster, normalized to the scale radius rg^ale- For a-U the 
simulated clusters rgcalc ~ 1 pc). Cross-hatched histogram (blue 
on the web): radial projected distribution of observed ULXs from 
the closest star-forming region (from the sample in table 9 of 
Berghea 2009). Hatched histogram (green on the web): radial pro- 
jected distance of observed X-ray sources (with X-ray luminosity 
Lx > 5 X 10'^® erg s~^) from the closest star cluster (considering 
M 82, NGC 1569, NGC 5253, from figure 3 of Kaaret et al. 2004b). 
Cross-hatched and hatched histograms refer to right-hand j/— axis 
(i.e. the number of observed X-ray sources, A^x) and to the top 
X— axis (the offset, i.e. the projected distance of the X-ray source 
from the closest star cluster and/or star forming region, in pc). 
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Figure 3. Distribution of three-dimensional recoil velocities 
(vrcc) of the ejected MSBHs. 
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Figure 5. Semi-major axis a (left-hand panel), period P (central panel) and eccentricity e (right-hand panel) from simulations of MSBH 
systems. In all the panels: the empty black histogram represents the initial conditions, the filled histogram (red on the web) shows the 
sub-sample of the ejected MSBHs at time t = 10 Myr, whereas the hatched histogram (blue on the web) shows the sub-sample of the 
non-ejected MSBHs at time t = 10 Myr. 




Figure 6. Fractional variation of the semi-major axis a (left-hand panel), of the period P (central panel) and of the eccentricity e 
(right-hand panel) from simulations of MSBH systems. The filled histogram (red on the web) and the hatched histogram (blue on the 
web) show the sub-sample of the ejected MSBHs and that of the non-ejected MSBHs at t = 10 Myr, respectively. 




[Myr] 



Figure 4. Distribution of ejection times of MSBHs from the clus- 
ter (toj). 



Table 2. End states of the simulated MSBHs that, at t 
Myr, are still bound to the cluster. 



10 



State 


Percentage'' 


Ionization 


7.2 


Binary 


21.4 


Triple 


71.4 


Exchange 


0.0 



''The percentage refers to the total of MSBHs that are still 
inside the parent cluster at £ = 10 Myr (i.e. 14 cases). 



4 EVOLUTION OF THE ORBITAL 
PARAMETERS 

Fig. 5 shows the main orbital parameters of the simulated 
MSBH binaries, i.e., from left to right, semi-major axis a, pe- 
riod P and eccentricity e. In all panels, the black empty his- 
togram represents the initial conditions, whereas the ejected 
MSBHs (at t = 10 Myr) and the MSBHs remaining in the 
cluster are represented by the filled histogram (red on the 
web) and by the hatched histogram (blue on the web), re- 
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Figure 7. Distribution of the fractional changes of the Roche 
lobe ( )■ Hatched histogram (blue on the web): runs where 
the MSBH is still inside the parent cluster at t = 10 Myr. Filled 
histogram (red on the web): runs where the MSBH is ejected 
before time t = 10 Myr. 



spectively. The initial conditions illustrate what we men- 
tioned in Section 2, i.e. the fact that we consider only hard 
binaries. A large fraction of the systems have periods typi- 
cal of RLO-HMXBs binaries (P < 10 days). After 10 Myr, 
there is no big difference between ejected systems and bina- 
ries that remain inside the parent cluster. The only notice- 
able difference concerns the final distribution of eccentrici- 
ties: after 10 Myr, the eccentricity of ejected binaries is, on 
average, larger than the eccentricity of systems that stay in- 
side the cluster. This agrees with the expectation that close 
encounters increase the eccentricity of a binary (see, e.g., 
Sigurdsson & Phinney 1993; Rasio & Heggie 1995; Heggie & 
Rasio 1996; Colpi, MapeUi & Possenti 2003). 

In order to highlight the differences between ejected 
and non-ejected systems, it is better to look at the frac- 
tional change of the orbital parameters shown in Fig. 6. 
From left to right, the fractional change of semi- major 
axis (Aa/a = [a(t = 10 Myr) - a(t = 0)] /a(t = 0)), pe- 
riod (AP/P = [Pit = 10 Myr) - Pit = 0)] /P(t = 0)) and 
eccentricity (Ae/e = [e{t = 10 Myr) - e(t = 0)] /e(t = 0)) 
are shown. The main feature of non-ejected systems is 
that the fractional changes of semi-major axis and of pe- 
riod are extremely small (|^| < 0.02 and |^| < 0.025). 
This indicates that non-ejected systems did not undergo a 
strong interaction. Instead, fractional changes of a and P 
for ejected systems are much larger (0.01 < |^| < 0.54 and 
0.01 < 1-^1 < 0.7), indicating a strong perturbation before 
the ejection. Furthermore, the semi-major axis decreases in 
most of the systems: ^ < for 91 and 85 per cent of ejected 
and non-ejected binaries, respectively. Again, this is consis- 
tent with the fact that the simulated binaries are hard and 
tend to become harder (Heggie 1975). It is interesting to 
note that only one ejected binary becomes wider and — is 
quite large (^ 0.5). However, this system is quite peculiar, 
as it was by far the softest already in the initial conditions 
(a(i = 0) = 2.8 A.U.), it has been ejected late in time (at 



t — 9.25 Myr), it is not far from the centre of the cluster 
(projected distance r = 6.8 pc at f = 10 Myr) and is one of 
the systems that underwent exchange of the companion. 

The right-hand panel of Fig. 6 shows the fractional 
change of eccentricity and strengthens the evidence that 
ejected binaries become more eccentric than non-ejected 
ones. 

In general, three-body encounters change significantly 
the orbital parameters of MSBH binaries: this may be cru- 
cial for their evolution into RLO-HMXB systems. In fact, 
three-body encounters may have the same effects as natal 
kicks: (i) they reduce the semi-major axis and can bring 
into RLO systems that, otherwise, would be too wide; (ii) 
they increase the eccentricity and may switch-on the RLO- 
HMXB phase, since a highly eccentric system is likely to be 
in RLO at periastron (and frictional forces rapidly circular- 
ize the orbits, leading to a stable RLO). 

Thus, dedicated runs including stellar and binary evo- 
lution are required, to quantify the effect of gravitational 
encounters on mass transfer and to determine the onset of 
a RLO-HMXB phase in MSBH systems. In this paper, we 
can only try to estimate the effects caused by changes in the 
orbital parameters, by calculating the Roche lobe (Pl). We 
adopt the most common approximation for the Roche lobe 
for circular orbits (see, e.g., Tauris & van den Heuvel 2006): 



a • 0.6g2/3 + ln(l-hqi/3)' ^ ^ 

where q — m2/mBH (where rriBH is the mass of the 
MSBH and m2 the mass of the companion). Fig 7 shows 
the fractional changes of the Roche lobe (APl/Pl = 
[PL(t = 10 Myr) - Ri^{t = 0)] /RL{t = 0)) for non-ejected 
(hatched histogram) and ejected MSBH binaries (filled his- 
togram). We note that the Roche lobe of MSBHs remaining 
inside the parent cluster is essentially unchanged, whereas 
the Roche lobe of ejected MSBH binaries may significantly 
decrease (|-^^| < 0.5). This is a direct consequence of the 
shrinking of the binary after a three-body encounter (see 
the left-hand panel of Fig. 6), combined with the possible 
change of m2/mBH through dynamical exchanges. Dynami- 
cal exchanges may be crucial for entering the RLO regime, 
as they generally lead to binary systems with more massive 
companion stars. This suggests that the orbital shrinking 
may switch-on RLO in ejected binaries, although our simu- 
lations cannot be conclusive, because of the small statistics 
and the absence of stellar evolution calculations^. 



5 COMPARISON WITH IMBHS 

For comparison with the MSBHs, we run 20 simulations 
of clusters having the same properties as before, but host- 

^ Because of the absence of stellar evolution in our simulations, 
we cannot estimate the radius of the companion star r, . There- 
fore, we do not know the value of the ratio i?L/r« and we cannot 
quantify whether a simulated system undergoes RLO. However, 
for reasonable assumptions about the stellar radius-mass relation 
(see, e.g., Prialnik 2000), Ri^/r* 0.8 — 20 in the initial condi- 
tions. Ri^/rtf remains almost constant for the retained binaries. 
The number of ejected binaries with Rj^/r, ^ 1 is four at t = 10 
Myr (only two of them had Rj^/rt ^ 1 at i = 0). 
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ing a 300 M0 IMBH. The adopted mass is quite high for 
some scenarios of IMBH formation. For example, Portegies 
Zwart & McMillan (2002) indicate that the mass of an IMBH 
formed via runaway collapse should be about 0.1 per cent 
of the mass of the parent cluster. Actually, the prediction 
by Portegies Zwart & McMiUan (2002) should likely be con- 
sidered an upper limit, as it does not include recipes for 
mass loss during the collapse (e.g., Gaburov, Lombardi & 
Portegies Zwart 2010). Thus, our assumption for the mass 
of the IMBH is somewhat unrealistic. On the other hand, 
it would be highly expensive to simulate a statistically sig- 
nificant sample of 3 x 10^ M© clusters with the required 
accuracy. Furthermore, the density in the core of the sim- 
ulated clusters (10^ — 10^ stars pc~^) is not far from the 
density of more massive clusters and the dynamical evolu- 
tion of the IMBH is expected to be quite similar in the two 
environments. 

Fig. 8 shows the two-dimensional projected distance of 
IMBHs from the centre of the parent cluster at t = (empty 
histogram) and t = 10 Myr (filled histogram). The bottom 
panel of Fig. 9 shows the projected distance of IMBHs at 
t = 10 Myr (hatched histogram), in comparison with the 
projected distance of MSBHs at t = 10 Myr (filled his- 
togram). First, we note that only one IMBH (5 per cent 
of the sample) has been ejected from the cluster after t = 10 
Myr. This is not surprising, as we expect recoil velocities 
for IMBHs to be at least a factor of 3 lower than those for 
MSBHs. We note that the distance of IMBHs from the cen- 
tre of the parent cluster at t = 10 Myr is larger than in 
the initial conditions. This was true also for MSBHs (see 
Fig. 2) and is likely due to the fact that three-body en- 
counters push into the outskirts of the cluster even those 
BHs that were not completely ejected (as predicted, e.g., by 
Kulkarni et al. 1993). We stress that this migration of BH bi- 
naries to the outskirts does not affect the global distribution 
of stars in the cluster, as the overall profile does not change 
significantly (see Fig. 1). It must be noted that MSBHs and 
IMBHs that are pushed into the periphery of the parent clus- 
ter might be more easily stripped from it, by external tidal 
fields (e.g. the one of the host galaxy). On the other hand, 
the escape of IMBHs and MSBHs from the parent cluster, 
because of tidal fields, implies (i) longer timescales than di- 
rect three-body ejection (likely longer than the lifetime of 
a massive companion star); (ii) the lack of recoil velocity. 
Thus, the simulations confirm that IMBHs can hardly be 
ejected from the parent cluster by three-body encounters. 



6 COMPARISON WITH DATA 

Fig. 2 compares the simulations of MSBHs with two differ- 
ent observational samples (see also Fig. 9, for a comparison 
of the two observational samples with the simulations of 
both MSBHs and IMBHs). The hatched histogram (green 
on the web) shows the offset of bright X-ray sources (with 
X-ray luminosity Lx > 5 x 10^^ erg s~^) from the closest 
star cluster in a sample of three nearby starburst galaxies: 
M82, NGC 1569 and NGC 5253 (from fig. 3 of Kaaret et al. 
2004b). The observed displacements reported by Kaaret et 
al. (2004b) range from ~ 10 pc to ~ 1000 pc, in good agree- 
ment with the projected positions of ejected MSBHs in our 
simulations (filled histogram). We note that displacements 




Figure 8. Empty black histogram: two-dimensional radial pro- 
jected distribution of IMBHs at time t = 0. Filled histogram 
(red on the web): two-dimensional radial projected distribution of 
IMBHs at time t = 10 Myr. j;— axis: number of simulated IMBHs 
(^imbh); i:— axis: two-dimensional distance of the IMBH from 
the centre of the parent cluster r, normalized to the scale radius 
^scalc- For all the simulated clusters rgcalc ~ 1 pc. 



smaller than 10 pc have been set to 10 pc by Kaaret et al. 
(2004b), because their spatial uncertainty is about 10 pc. 
Therefore, nothing can be said about X-ray sources inside 
the parent cluster, on the basis of Kaaret et al. (2004b). 

The cross-hatched histogram (blue on the web) shows 
the offsets of a sample of ULXs from the closest region of 
star formation (data come from table 9 of Berghea 2009). 
Although the offsets reported by Berghea (2009) refer gener- 
ically to star forming regions (which cannot be distinguished 
from YSCs), the superposition between the position of simu- 
lated ejected MSBHs and the data is remarkable. Both peak 
at ~ 100 pc and range from ~ 4 pc up to ~ 1000 pc. Only the 
position of a single ULX (IC 342 X-3) in Berghea (2009) is 
found to be completely coincident with the centre of a star 
forming region. This might be at odds with the fact that 
more than half of the simulated MSBHs at t = 10 Myr are 
still inside the parent cluster, but there are various possible 
explanations for this difference. 

First, data by Berghea (2009) represent the offset be- 
tween the ULX and the centre of the closest star forming 
region, identified on the basis of the colour; but each star 
forming region may host more than one YSC and/or OB 
association. Thus, the position of the ULX might be coinci- 
dent with that of a YSC enclosed in the star forming region, 
but not identified by Berghea (2009). Berghea (2009) high- 
lights that about half of the ULXs in his sample are inside 
the star forming region, although displaced with respect to 
its centre. Furthermore, table 8 of Berghea (2009) indicates 
that, out of a sample of 31 ULXs, 6 have no detected coun- 
terpart (2 of them do not even have a star forming region 
within 1 kpc), 12 have single-star candidates as counterparts 
and 13 have counterparts that are brighter than a single 
05V star. Among the latter, 5 ULXs are associated with 
ultra-violet (UV) emission consistent with > 14 05V stars. 
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7 CONCLUSIONS 




0.01 0.1 1 10 100 1000 
offset [pc] 



Figure 9. Bottom panel: two-dimensional radial projected dis- 
tribution of IMBHs (hatched histogram) and of MSBHs (filled 
histogram, red on the web) at time t = 10 Myr. Bottom y— axis: 
number of simulated BHs (A^bh); a;— axis: two-dimensional dis- 
tance of the BH from the centre of the parent cluster (offset) 
in pc. Top panel: the cross-hatched histogram (blue on the web) 
shows the radial projected distribution of observed ULXs from 
the closest star-forming region (from the sample in table 9 of 
Berghea 2009); the hatched histogram (green on the web) shows 
the radial projected distance of observed X-ray sources (with X- 
ray luminosity Lx > 5x 10^^ erg s~^) from the closest star cluster 
(considering M 82, NGC 1569, NGC 5253, from figure 3 of Kaaret 
et al. 2004b). Top j/— axis: the number of observed X-ray sources 
(Nx); X— axis: projected distance of the X-ray source from the 
closest star cluster and/or star forming region, in pc. 



which possibly means YSCs similar to those we simulated 
(although the UV excess may be due in part to contamina- 
tion from the accretion disc). In particular, Berghea (2009) 
highlights that Holmberg IX X-1 is located at the edge of a 
stellar cluster with estimated mass of 10^ Mq (Grise, PakuU 
& Motch 2006), NGC 1313 X-2 is located in the periphery 
of a relatively young cluster (Zampieri et al. 2004), IC 342 
X-3 is inside the nuclear cluster (Boker, van der Marel & 
Vacca 1999) and NGC 4449 X-7 is associated with a com- 
pact young cluster (Gelatt, Hunter & Gallagher 2001; Soria 
et al. 2005). Therefore, the data by Berghea (2009) do not 
necessarily exclude the association of a fraction of ULXs 
with YSCs. 

On the other hand, there may be intrinsic reasons why 
HMXBs from MSBHs are suppressed inside star clusters. 
For example, a MSBH binary that remains inside the clus- 
ter continuously undergoes three-body interactions (see Ta- 
ble 2) and cannot evolve unperturbed. This might strongly 
affect the mass transfer (see, e.g., Mapelli et al. 2004, 2006 
for the effect of dynamical perturbations onto blue straggler 
formation from mass-transfer binaries). These effects might 
be very important and need to be studied with dedicated 
binary-evolution simulations. 



In this paper, we studied the dynamical evolution of hard bi- 
naries hosting MSBHs with massive companions. We showed 
that a large fraction (~ 44 per cent) of MSBHs is ejected 
from the cluster on a short timescale (< 10 Myr), because 
of close interactions. This result is very different from what 
the simulations indicate for IMBHs; only a small fraction of 
IMBHs (~ 5 per cent) is completely ejected from the parent 
cluster. 

All the ejected MSBHs retain their companion, even af- 
ter the ejection. Furthermore, at t ~ 10 Myr, the offset of 
the ejected MSBHs with respect to the centre of the par- 
ent cluster is consistent with the observations of Kaaret et 
al. (2004b) for HMXBs and with those of Berghea (2009) 
for ULXs. These results also suggest that dynamical recoil 
for MSBHs has similar effects to natal kick for stellar BHs. 
Thus, assuming that MSBHs born from direct collapse have 
no (or negligible) natal kick, the dynamical recoil is the only 
efficient mechanism to eject MSBH binaries from their par- 
ent cluster. 

In some cases (especially for IMBHs), the recoil due to 
three-body encounters is not sufficient to promptly eject the 
BH, but the BH is pushed into the outskirts of the clus- 
ter. If this occurs, the BH might escape at a later stage, 
because of the influence of the galactic tidal field or cluster 
evaporation /disruption . 

Also, three-body interactions might affect the orbital 
parameters of the MSBH binaries. The semi-major axis and, 
consequently, the period change sensibly after the ejection. 
On average, the binary hardens, and experiences a frac- 
tional decrement of the semi-major axis of 0.01 — 0.5. Ejected 
MSBH binaries tend to become more eccentric. These effects 
on the orbital parameters of MSBH binaries might have im- 
portant consequences on the occurrence of mass transfer. In 
fact, a decrease of the semi- major axis might bring into RLO 
systems that were, otherwise, too far apart. 

Recent studies (Mapelli et al. 2009; Zampieri & Roberts 
2009; MapeUi et al. 2010a,b) propose that MSBH binaries 
might power a fraction of bright HMXBs and ULXs. On the 
other hand, LIO indicate that MSBH binaries can hardly 
become RLO-HMXBs, because of the lack of natal kick 
in their formation pathway. In this paper, we suggest that 
dynamical interactions may have similar effects as the na- 
tal kick. This argument is strengthened by the agreement 
between the radial distribution of simulated MSBHs and 
the offset observed in many ULXs and/or HMXBs with re- 
spect to nearby star clusters and/or star forming regions. We 
also suggest that dynamical interactions can change the or- 
bital parameters of MSBH binaries, favouring mass transfer. 
Thus, it will be crucial to assess whether dynamical inter- 
actions may effectively lead MSBH systems to pass through 
RLO-HMXB phase. For this reason, dedicated binary evo- 
lution calculations, coupled with dynamical simulations, are 
definitely needed. 

Furthermore, the aforementioned papers (Mapelli et al. 
2009; Zampieri & Roberts 2009; LIO; MapeUi et al. 2010a,b) 
consider only the formation of massive BHs as remnants of 
metal-poor stars. However, there might be different path- 
ways for the formation of massive BHs, even independent 
of the metallicity. For example, in dense star clusters, mas- 
sive BHs are expected to form through (single or multiple) 



Dynamics of MSBHs in YSCs and ULXs 9 



mergers of hard binaries, as a consequence of stellar evolu- 
tion (e.g., Belczynski et al. 2004), 3-body encounters (e.g.. 
Miller & Hamilton 2002) and/or (in case of binaries of com- 
pact objects) gravitational wave emission. In addition, the 
runaway collapse process (e.g., Portcgics Zwart & McMil- 
lan 2002; Giirkan, Freitag & Rasio 2004; Freitag, Giirkan & 
Rasio 2006), which is expected to lead to the formation of 
IMBHs, can be interrupted when the BH mass is still < 100 
Mq, because of the dynamical ejection of the BH seed (e.g., 
Portcgics Zwart & McMillan 2000; Miller & Colbert 2004) 
or of strong mass losses in the merger phase (e.g., Gaburov, 
Lombardi & Portegies Zwart 2010). Such massive BHs born 
from an interrupted runaway collapse or from binary merg- 
ers need further investigation. 
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